Time evolution of Wigner function in laser process derived by entangled state 

representation 



o 
o 

(N 



^ ■ 

^— > ■ 

c ■ 

s : 



> 
o 
o 

(N 

cn 
o 
o> 
o 



X 



Li-yun Hu 1,2 *and Hong-yi Fan 2 
College of Physics & Communication Electronics, 
Jiangxi Normal University, Nanchang 330022, China 
2 Department of Physics, Shanghai Jiao Tong University, Shanghai, 200030, P.R. China and 
* Corresponding author. E-mail addresses: hlyun2008@126.com 

Evaluating the Wigner function of quantum states in the entangled state representation is intro- 
duced, based on which we present a new approach for deriving time evolution formula of Wigner 
function in laser process. Application of this fomula to calculating time evolution of photon number 
is also presented, as an example, the case when the initial state is photon-added coherent state is 
discussed. 



I. INTRODUCTION 

One of the major topics in Quantum Statistical Me- 
chanics is the evolution of pure states into mixed states 
P, 0] ■ Such evolution usually happens when a system is 
immersed in a thermal environment, or a signal (a quan- 
tum state) passes through a quantum channel, and is 
described by a master equation. Alternately, descrip- 
tion of evolution of density matrices p can be replaced 
by its Wigner function's evolution in phase space [1, 0| . 
The partial negativity of Wigner function can be consid- 
ered as an indicator of nonclassicality of quantum state. 
On the basis of the entangled state representation and 
the thermo field dynamics we present a new approach 
for deriving time evolution formula of Wigner function 
in amplitude-damping channel and laser process. Ap- 
plication of this fomula to calculating time evolution of 
photon number is also presented, as an example, the case 
when the initial state is photon-added coherent state is 
discussed. 



II. WIGNER FUNCTION FORMULA IN 
THERMO ENTANGLED STATE 
REPRESENTATION 



We begin with briefly reviewing the thermo entangled 
state representation (TESR) . On the basis of Umezawa- 
Takahash thermo field dynamcs (TFD) @, i, we have 
constructed the TESR in doubled Fock space [1, Q , 
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equations of \rf), 

(a - a f ) \r)) = Tj |??) , (a f - a) \rj) = rf \rf) , 

fa| (at - a) = rf ( V \ , (r,\{a-^)= V (v\- (3) 

Note that Ua — a'), (<r — 5)1 = 0, thus \rf) is the 
common eigenvector of (a — aJ) and (5 — a'). Us- 
ing the normally ordered form of vacuum projector 
|0, 5) (0, 6| = : exp (— a^a — a^a) : , and the technique of 
integration within an ordered product (IWOP) of opera- 
tors [13, [IH 0] > w e can easily prove that \rj) is complete 
and orthonormal, 
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It is easily seen that \rj = 0) has the properties 
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Note that density operators p(a> ,a) are defined in the 
real space which are commutative with operators (a* ,a) 
in the tilde space. 

Next, we shall derive a new expression of Wigner func- 
tion in the TESR. According to the definition of Wigner 
function [l3|, [3] of density operator p. 



\r)) = D(ri)\r, = 0), D( V ) = 



(2) 



where D (77) is the displacement operator, a) is a fictitious 
mode accompanying the real photon creation operator a) , 
|0,6) = |0) 1 6) , and |5) is annihilated by a, [a, a^] = 1. 
Operating a and a on \rj) in Eq.((T]) we obtain the eigen- 



W (a) = Tr [A (a) p] . 



(7) 



where A (a) is the single-mode Wigner operator [13j |. 
whose explicit form is 



A (a) = -: r 

7T 



-2(a + -Q*)(a-a) . 
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By using (n| ih) = S n-m and introducing Ip) = p\I) we 
can reform Eq.© as 

oo 

W (a) = (n, h\ A (a) p \m, in) 
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(t7 = 0| D (2a) (~l) a ' a \p) 
(rj = -2a\(-ir ia \p) 
(€ - 2a| p) , 



where |£) is defined as 

l0e=, = (-i) o, °l»? = -0 



exp(--|C| 2 +ea t +rS t -a t 5 t ) |0,6) 
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= D(£)e-° a |0,6) 
It can be proved that 
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\z. z*) , we have 
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W\ z )( z \ (a) = - (0, 0| exp (— 2|«| 2 + 2a* a) + 2aa - aa) \z, z* 
1 
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-2\a- z\ 



(16) 



Further, using Eq. (lll|) and the completeness of (r)\ in 
Eq.©, we can reform Eq.© as 

W («)=/$ = 2a \ M = / e "*^ Q "* M <°> 

(17) 

Once (t)\ p) is known, one can calculate the Wigner func- 
tion by taking the Fourier transform of (r?| p). Eqs. © 
and (flU|) are two ways accessing to Wigner function, we 
can use either one to derive Wigner functions. 



III. EVOLUTION FORMULA OF WIGNER 
FUNCTION FOR AMPLITUDE DAMPING 
CHANNEL 



(m = 2 ex P I — 2 — 



(ii) 



a Fourier transformation kernel, so |£) can be considered 
the conjugate state of \rj) , which also possess orthonor- 
mal and complete properties 

(£1 = 1, {?\t)=*6tf-Q5(f-?). 

7T 

(12) 

Eq.© is just a new formula for evaluating the Wigner 
function of quantum states: by calculating the overlap 
between two "pure states" in enlarged Fock space rather 
than using the ensemble average in real mode space. 

For example, for number state \n) (n\ , noticing 
\ri) (n\ I) — \n,n), and the generating function of two- 
variable Hermite polynomial [H, HBl H mn (x,y), 

V ——IT (a?, y) = exp [-i£ + tx + t'y] , (13) 



we see 

W ln){nl (a) = - (n,n\ £ =2q ) = -j-e-^X.n (£,£* ) 
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in the last step in Eq. (fT4")) we have used the relation be- 
tween H m n (x,y) and Laguerre polynomial L m (x) [l7| . 



(-1)" 

L n (xy) = 7—H niTl (x,y) 



(15) 



In this section, we consider Wigner function's time evo- 
lution in the amplitude decay channel (dissipation in a 
lossy cavity) described by the following master equation 

M 



= k (2apat — at ap — pat a) , (18) 

where n is the rate of decay. In Ref. [2l[ we have re- 
formed (fT8l) as 



(19) 



— \p) = K (2aa — at a — at a) \p) , 
thus the formal solution of Eq. (fT5]) is 

\p(t)) = e Kt ( Qa - atQt + l) e (l-e 2Kt )(a t -a)(a-5 t )/2 

(20) 

Then projecting Eq. (f20|) on (77I, and noticing 
exp hit (aa — aW)] being the two-mode squeezing 
operator, 



Pa) 



(77 1 exp hzt (aa — at at 



)] = e- Kt (rje-^l 



as well as Eq.©, we obtain 

- p -hTH 2 



(r,\p(t)) = e-^ ( V e- Kt \ Po ), 



(21) 



(22) 



where T = 1 - e~ 2Kt . Substituting Eq.([22D into Eq. (|T7|> . 
we derive the Wigner function at time t 



W(a,t) 



- [ ^JL„a*,,-a n --^T\7 1 



27T 2 



' 2 {v e ~ Kt \ Po) ■ (23) 



Similarly, for coherent state \z) (z\ (\z) = exp(— \z\ /2 + 
zat) |0)) [H, due to \z) (z\ I) = D (z) D (z*) |00) = 



Inserting the completeness relation (fT2")) into Eq. (f2"3")) and 
noticing Eqs.© as well as (fTTj) , we can reform Eq. 
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W (a, t) 



IV. EVOLUTION FORMULA OF WIGNER 
FUNCTION FOR LASER PROCESS 



(Ve- Kt \0 Po) 

IT J 27T Z 

r j2 oj2 We now generalize the master equation to the case of 

/ 1 e -% h\ 2 +v{a'-i3*e- rit )+7 1 '-(i3e- K - t -a) w ^I^ser theory. The mechanism of laser is described by the 

J ft following master eauation 



d 2 (3 



• exp 



T 



\a — j3e 



W(J3,0) 



(24) 



where W (/?, 0) is the Wigner function at initial time, and 
we have used the following integral formula [ItJ 



^exp ((\z\ 2 +(z + V z*^ =~e-^,Re(C)<0. 

(25) 

Eg. (f2"4"| is the expression of time evolution of Wigner 
function for amplitude damping channel. 

For example, for the photon-added coherent state 
C m a) m \z), where C m = [m\L m (— |z| 2 )] -1 is the nor- 
malization factor, the initial Wigner function W (/?, 0) 
is given by [22[ 

W(J3,0)= [ ' L m {\2f3-z\ 2 ). (26) 

7rL m (- \z\ ) 



Substituting Eq.d26j) into Eq. ([24)l and using Eq.([T5j) as 
well as the another generating function of H m ^ n (x, y) , 



H m ,n (x,y) 







im+n 



g T mg T tn 



exp [—tt + tx + t y] 



we have 
W(a,t) 



t=t'=0 

(27) 



2 e -2(\z\ 2 + ±\ a \ 2 



d 



2m 



-TT —TZ — Z T 



T nm\L m (-\z\ 2 )dr m d^ 



d 2 p 



7T 



■ exp 



2 
T 



+ 2(3[z* + —e- 



-2/3* ( 



o 



z H e 

T 

— 2\a. — ze 



— Mtia 



T 

Q2m 



t=t'=0 



exp [(1 - 2e- 2tK ) tt' 



nm\L m (- \z\ 2 ) dT m dr' m 
+ [(1 - 2e- 2tK ) z + 2ae~ tK ] r 
+ [(l-2e~ 2t *)z* +2a*e- tfi ]r'] r=T , =0 . (28) 

With use of a scaled transformation in the right-hand 
part of Eq. ([25)1 we finally get 



W(a,t) 



(l-2e- 2Kt ) m 
wL m (- \z\ ) 



x L, 



\2ae~ 



(l - 2e- 2Kt 



I 2 T 



1 - 2e 



-2ret 



(29) 



which is the analytical expression of the time evolution 
of Wigner function for any number (m) photon-added 
coherent state in photon loss channel [23j. In particular, 
when t = 0, Eq. ([2"9"[) just reduce to Eq. 



following master equation 

^rr - = g [2a 1 p (t) a — aa) p (t) — p (t) oa'] 
+ k [2ap (t) a 1 — a) ap (t) — p (t) a* ail 



(30) 



where g and k are the cavity gain and the loss, respec- 
tively. Eq. (|3"D|) reduces to Eq. tfTS]) when 5 = 0; while for 
g — ► nn and K — * k (n + 1) , Eq. (|30p becomes 



= k (n + 1) (2apa) — a) ap — pa) a) 



(2a) pa — aa) p — paa' ) 



(31) 



which corresponds to the master equation in thermal en- 
vironment [20l | . 

Similar to the way of deriving Eq. (|22[) . we have derived 
in Ref. [H[ 



\p (t)) = exp \(aa — a'a) + l) (/t — g) t 
"(/c + s)(l-e 2 C"-»>*) 



x exp 



2 (k - 5 ) 



(a^ — a) (a — 



Thus the matrix element (77 1 p (i)) is given by 
A, 



(r)\ p(t)} = exp 



where 



.4 



K + g 
k- g 



1 



rje 



-{K-g)t 



Po) 



= _2( K -g)t 



Po) 
(32) 



(33) 



(34) 



According to Eq. ([13]) the Wigner function's evolution for 
Laser process is given by 



W(a,t) 
d 2 rj 



2tt 2 

d 2 id 2 r, 
d 2 £d 2 rj 



e -2-\v\ 2 +a*V-arf ^g-(K-ff)* 



Po) 



e --\ v \ +a * v - av * /-( K - g) t w ^ 



|.| ?) |= +7) ( a *_ | 3* e -C«- 9 )t) + ^*( / 3e-(~-»)'-a) w /^ ^ 



d 2 /3 



■ exp 



a-(3e- {K - 9)t 



W(0,O), (35) 



where we have used Eq. (|25l) . In particular, when g 
0, Eq.([3"5"[) reduces to Eq. (IM[) . For g — » ren and k - 
k (n + 1), leading to ^4 = (2n + 1) T, Eq.([55[l becomes 



W(a,t) 



(2n+l)T J tt 



J20 c-/3e-" 
^W{P,0)e- 2 C2« + 1)T 



(36) 
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W (a, t) = 2e 2Kt j d 2 (3W T ((3) W (e K *(a - Vf/3), o) , 

(37) 

2|g| 2 

where Wt (/?) = ^a+i) e 2fi+1 i s the Wigner function of 
the thermal state with mean photon number n. 

Similar to the way of deriving Eq. (|29|) . when the initial 
state is C m a^ m \z) , substituting Eq. (l2l))) into Eq. (f3"5)) we 
have 



o~C-2\0\ 2 



irL m {~\z\ 2 ) (2nT + l) 



(nT + l) m 



A 



(38) 



where 
A= 1 

B = 
C = 



-iKt 



IT 



(2Tn + 1) (n + 1)' 



ftT+l 



-0 



nT + 1 " ' (2fiT+l) v /(n + l)T 
I — -\ze Kt \ + 4T 2 n 2 ' ' 



2nT + 1 \Tn+l 
2e~ Kt (Tn + 1 



2nT+l 



(Z0* +(3z*) 



(39) 



)|2 



In particular, when n = 0, leading to A = T 
-± T ((l-2e- 2Kt ) z* + 2e~ Kt f3*) , and -C - 2 |. , 

-2 |/3 - ze~ tK | 2 , thus Eq.® reduces to Eq.flU]). 

Eq. (f3"5|) manifestly shows that the Wigner function 
of C m a im \z) in thermal environment is closely related 
to the Laguerre polynomials. In addition, due to 
L m (^— \x\ 2 ^j > 0, so C m > 0, thus it is easily seen that 
when A > 0, which means the condition 



1, 2(n+l) 
nt ^ nt r = — In . 

2 2n + 1 ' 



(40) 



the Wigner function (l38|) is always positive-definite. 
Thus we emphasize that for any values of m, when the 
condition (|40| is satisfied, the Wigner function has no 
chance to be negative. 



V. TIME EVOLUTION OF PHOTON NUMBER 
FOR THE LASER PROCESS 

Next we consider the photon number (PN) of density 
operator p for the laser process. According to the TFD, 
we can reform the PN p (n) = (n\ p \ n) as 

oo 

P ( n ) = ( n \ P W = y~] (n, n\ p \m, in) 



(n,n\p\I) = (n,n\ p) 



(41) 



thus the PN is converted to the matrix element (n,n\ p) 
in thermo dynamics frame. Then using the completeness 
of (£| and Eq.Q as well as Eq. (fl"4|) . we see 

P(n)= I ^ (n,n\ ® (£\ p) 



d 2 ti (n,n\0W(a = £/2) 



= 4tt / dW|„ )W (a)ff(a), (42) 

one can see this formula also in [l|, [13]. Thus one can 
calculate the PN by combining Eq.J35|) and (|l2l . 

Now we evaluate the PN of the above decoherence 
model in Eq.JH). Substituting Eq.([35|) into Eq.(g2]), wc 

see 



p(n) = j I d 2 0W(0,O)G(p) 



(43) 



where 



G(fi) = j d 2 aW ln){nl ( a)e -ih-^- ( "- 9) T-^l 
d 2 a 



(-!)» a -^L n (A\a\ 2 )e-M 



(44) 



Using Eqs.([25|) and (J2TJ) we can evaluate Eq.(g4]) as 



x exp 



- 2 ^±l\a\ 2 + 2a(T+ P * 



A 



-2a* t' 



t=t'=0 



^4 e Ji+i iPl Qn+n 

2(A+l)n\ dr n dT> 



■ exp 



I- A 
1 + A' 



23*6-^-^ , 2l3e-^-^ 1 
H : : T + T 



A + l 



A + l 



■ (45) 

t=t'=Q 

After making some scaled transformations, wc finally ob- 
tain 



2{l + A) n+1 n\ 



d n+n -r'r+ 

x e Vi-a 2 

Q T nQ T ln 



e ±j r . 



t=t'=0 

f 4\0\ 2 e~ 2 ^9)t s 



2 (1 + A) 



n+l 



I- A 2 



(46) 



(A + 1) 

' L "< i_ A 2 l/91 a K09,O), (47) 



Substituting Eq.gHD into Eq.(g31) yields 
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which is a new formula for calculating the photon number 
distribution of the open system in enviornment. From 
Eq. (j47| it is easily seen that once the Wigner function of 
initial state is known, one can obtain its photon number 
distribution by performing the integration in Eq.(|47|). 

In particular, when g = 0, A = 1 - e~ 2Kt = T, Eq. ([47|) 
reduces to 



p(n) 



4(-l)"e 



(2e 2Kt - 1) 



Tt + l 



d 2 (3e 



rl/3| 2 



x L r , 



4e z 



2e 2 



1 



(48) 



which corresponds to the photon number of density op- 
erator in the amplitude damping quantum channel. 

While for g — > nn and k — ► «; (n + 1), , Eq. (|4"T)) be- 
comes to 



p(n) 



A {A 



{A 

x L n 





l) n 


1) 




4e 


-2ret 


1 - 


-.4 2 



(49) 



where .4 = (2n+l)T = (2n + 1) (l - e- 2Kt ) . Eq.® 
corresponds to the photon number of system interacting 
with thermal bath. 

For example, we still consider the photon-added co- 
herent state field. Substituting Eq. (|2l)|) into Eq. flTj) and 
uisng Eqs.([25]) and (gT]) yields 

p(») = Ne~^ j ^L m (W-z\ 2 )L n |/?| 2 



x exp 



2{z(i* + (3z*) -2 1 + 



Ne -2\z\ 2 



d 



2m 



,-2( K -g)t 

A + l 

Q2n 



-w —z V 



Q v mQ v lm Q T nQ T l 



X e 



where we have set 

2 - (i 



A = — , er 
A 4 



e -(«-ff)t 
v 7 ! - A 2 "' 



(51) 



and 



(A+l)" +1 k| 2 ) 



-2{ K -g)t 

A+l • 



(52) 



Further expanding the exponential 
exp [cijui/ + (Act — 1) tt'] , we finally obtain 

, , iVA 2 "e^l z l 2 £5 m!n! [ W (Aa-l)/A 2 l 

p(n) = > 

1ii{-uS) n - m ^ o /!fc![(m-0!(n-fc)!] 

2 



item 



x H, 



-Ln-k (ty/tJZ,ly/UZ 



(53) 



In particular, when g = 0, leading to .A = a; = T, a 
^ x l T i , /x = 2- e -^t , Act = 1, and A = y/l + T, thus 



p(n) 



ml (1 - w)" A uj m - n e- e ~ 



*i\ L m {~\z\ 2 )^ ll[(m-m 

|2 



X |i? m -/,„ (iy/wz, iy/wz*) I 



(54) 



■ exp 



-2/i |/3| + 2 (<r + 2* + v) (3 



+2(aT' + z + v')f3*} v=v , =T = T , =0 
_ N {-l) m+n j_^\ z? d 2m 



2-2fi 

-e n 



& 



2n 



2\im\n\ dv m dv' m dT n dr' n 

x exp [ojvv' + (Act — 1) tt' + A (tv' + vt') 
+c (*V + zv) + A (zr + «V)] w=r=T , =0 , (50) 



which concides with Eq.(43) with idea detection effi- 
ciency in Ref. [23] . 

In sum, by virtue of the thermo entangled state rep- 
resentation that has a fictitious mode as a counterpart 
mode of the system mode, we have derived the relation 
between the Wigner functions at t time and the initial 
time when quantum system interacts with envoirnment, 
such as decoherence, damping and amplification. As an- 
other quantity describing quantum system, the formula 
of photon number distribution has also been derived, 
which can be evaluated by performing an integration for 
the initial Wigner function. Our deriviations seem more 
concise. 
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